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The escape of a small current-biased Josephson tunnel junction from the zero voltage state in 
the presence of weak microwave radiation is investigated experimentally at low temperatures. The 
measurements of the junction switching current distribution indicate the macroscopic quantum tun- 
neling of the phase below a cross-over temperature of T* ~ 280 mK. At temperatures below T* 
we observe both single-photon and multi-photon transitions between the junction energy levels by 
applying microwave radiation in the frequency range between 10 GHz and 38 GHz to the junction. 
These observations reflect the anharmonicity of the junction potential containing only a small num- 
ber of levels. 



At low temperatures and small damping the dynamics 
of a current-biased Josephson junction is governed by the 
macroscopic quantum mechanics of the superconducting 
phase-difference across the junction (see e.g. Refs. [Q, ||] 
and references therein). Macroscopic quantum tunneling 
of the phase M , energy level quantization j|, |^] and the 
effect of dissipation [|| have been studied in detail 
Josephson junctions are solid-state quantum devices fab- 
ricated with integrated circuit technology. Their parame- 
ters can be adjusted in a wide range and can be well con- 
trolled. Josephson junction circuits have been proposed 
§, @, | and recently successfully tested §, El || |l| 
as qubits in quantum information processing [14|. 

In this letter we present experimental evidence of 
multi-photon transitions between the ground and the first 
excited state in a current-biased Josephson junction. The 
experiments have been performed below the cross-over 
temperature [jl5|, where the escape of the junction from 
a metastable state is dominated by quantum tunneling 
from the quantized energy levels. Using a high resolution 
measurement []l6| of the switching current JTtJ , we detect 
the multi-photon absorption by monitoring the decay of 
the junction from the zero-voltage to the finite voltage 
state. 

In the Stewart-McCumber model 19 1, the current 
biased small Josephson junction is modeled as a parti- 
cle of mass rn^ moving in an external washboard poten- 
tial = —Ej(j(p + cos</>), see Fig. 0, according to the 
equation of motion m^cj) + m < j ) (RC)~ 1 <f> + dU^/dcj) = 0. 
Here, the phase difference 4> across the junction repre- 
sents the position of the particle. The particle mass is 
given by = C(&o/2tt) 2 , where C is the junction ca- 
pacitance and $o is the flux quantum. Ej = ^ n I c /2n is 
the Josephson coupling energy with the critical current 
of the junction I c determining the depth of the poten- 
tial. The applied bias current / normalized as 7 = I / I c 
determines the tilt of the potential and the junction re- 
sistance R causes the damping proportional to the coef- 
ficient 1/RC. 

In the absence of thermal or quantum fluctuations and 
for 7 < 1, the junction is in the zero voltage state, cor- 
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FIG. 1: a) The Josephson junction energy U^(4>) calculated 
for the experimental parameters at the bias 7 = 0.995 (solid 
line). Numerically calculated energy levels (dotted lines) and 
the squared wave functions (dashed lines) are shown, b) 
Multi-photon transitions between the ground state and the 
first excited state. 



responding to the particle being localized in the the po- 
tential well. At finite temperatures T > 0, the parti- 
cle may escape from the well at bias currents 7 < 1 
by thermally activated processes O, ^(| or by quan- 
tum tunneling through the barrier |^ The rate at which 
both processes occur depends on the barrier height Uq = 
2Ej Ul-7 2 - 7 arccos( 7 )j w Ej 4^2/3 (1 - 7) 3 / 2 , the 
oscillation frequency of the particle at the bottom of the 
well 4 = ^/U"^{4> Q )/m^ = u p (1 - 7 2 ) 1/4 , see Fig. | 
and the damping in the junction. Here uj p = y/2nI c /QoC 
is the Josephson plasma frequency. At temperatures be- 
low the cross-over temperature T* ]l5f the quantum tun- 
neling rate dominates the thermal activation rate. The 
quantization of the energy of oscillations of the phase at 
the bottom of the well, see Fig. |l|a, has been observed 
both below [g and above T* J§. 

The experiments presented here were performed using 
a high quality 5x5 /1m 2 tunnel junction fabricated on an 
oxidized silicon waver using a standard Nb/Al-AlOz/Nb 
trilayer process. The junction had a critical current den- 
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sity of j c ph l.lkA/cm 2 , a capacitance of C ~ 1.6 pF 
and a subgap resistance of R > 500 57 at T < 2.0 K. For 
these sample parameters the expected energy level sep- 
aration is larger than 100 GHz at zero bias. The level 
width is small relative to the level spacing, due to the 
small damping. The predicted cross-over temperature 
T* is larger than 250 inK. 

The sample was mounted in an rf-tight sample box 
on the cold finger of a dilution refrigerator. The dc-bias 
leads were filtered with 7r-typc fccdthrough filters at room 
temperature, RC-filters at the 1 K-pot of the refrigerator 
and thermocoax filters at the sample box in order to re- 
duce external electromagnetic interference. A microwave 
signal was fed into the sample cell via a superconducting 
semi-rigid coaxial cable. To reduce the relative level of 
noise in the microwave signal, several stages of cold atten- 
uators of a total of —40 dB were used. We have verified 
that the power of all harmonics and subharmonics was 
at least 80 dB below the fundamental frequency power. 
For the switching current measurements [ p7| , the current 
was ramped up at a constant rate of / = 0.245 A/s with 
a repetition rate of 500 Hz. The switching current was 
determined by a measurement of the time delay between 
the zero-crossing of the bias current and the appearance 
of a voltage across the junction pl| . 

The switching current distribution P(I) of the sample 
in the absence of microwaves was measured in the tem- 
perature range between 4.2 K and 25 mK. These mea- 
surements indicate the thermal activation of the phase 
at high temperatures followed by a cross-over to quan- 
tum tunneling around T* ss 280 mK . In this letter 
we present measurements performed at T < T*. 

Microwaves in the frequency range between 10 GHz 
and 38 GHz were applied to the sample. While mon- 
itoring the P(I) distributions of the junction, the mi- 
crowave power P^ w was swept from low values, at which 
the P(I) distribution is not changed by the microwaves, 
to higher values for each chosen frequency. At negligibly 
small microwave powers the P(I) distribution is essen- 
tially determined by the unperturbed quantum tunnel- 
ing of the phase from the ground state of the well. If the 
microwave power is increased to substantially populate 
the excited level, the P{I) distribution becomes double- 
peaked. This double-peak structure smoothly varies with 
Pfj, w , as shown for v — 36.554 GHz in Fig. ||a. Further 
increasing the power, only the pronounced resonant peak 
is visible in the distribution. At this level of power the 
populations of the ground and the first excited state are 
equal, but the tunneling rate from the excited state is ex- 
ponentially larger than that from the ground state. Thus, 
the Pit) distribution is dominated by the resonant peak 
due to tunneling from the first excited state and the ini- 
tial peak in the distribution disappears, see Fig. | Due 
to the resonance excitation of transitions between the two 
levels, the switching current distribution at this level of 
power is more narrow than in the absence of microwaves. 



This fact proves that the measured P(I) distribution in 
absence of microwaves is not limited by noise in our ex- 
perimental setup and that below T* the escape indeed 
occurs due to quantum tunneling through the barrier. 

The bias current at which the resonant peak in the 
P(I) distribution appears depends strongly on the mi- 
crowave frequency. Most strikingly, we observe resonant 
peaks at similar or the same bias current for very different 
microwave frequencies. In Figs. |^b and c, two representa- 
tive density plots of the switching current distributions 
versus the applied microwave power are shown for the 
microwave frequencies 36.554 GHz and 18.399 GHz. For 
both frequencies the resonant peaks appear at almost 
identical bias currents. Both sets of data show the pro- 
nounced narrowing of the distribution at the resonance. 

The resonant bias currents I r , defined as the current 
at which the escape of the phase is maximally enhanced 
by the microwaves as indicated in Fig. || and Fig. [|, were 
extracted for all measured microwave frequencies v, see 
Fig. ||. It is clearly observed that the resonances fall into 
different groups as indicated by the dashed lines. We 
find that ratios of the resonance frequencies v for a fixed 
current I r are given with high accuracy by ratios m/n 
of two small integer numbers n and m, suggesting that 
the observed effect is related to multi-photon transitions 
between energy levels of the phase. 

In parabolic approximation of the potential, one ex- 
pects the energy level separation in a small Joseph- 
son junction to scale with the applied bias current as 
AE = hw p (l - (7 r // c ) 2 ) 1/4 - Therefore resonances with 
the external applied microwaves are expected to appear 
for nv = AE, where n is the number of photons ab- 
sorbed in the transition between two energy levels. Such 
multi-photon transitions between neighboring energy lev- 
els are quantum mechanically allowed [ plj in the anhar- 
monic potential for the phase of a Josephson junction 
due to the large diagonal matrix elements of the excited 
states. In Fig. || all data are fitted to the single for- 
mula (l/n)i/„(l - (Ir/Ic) 2 ) 1 ^, with v v = 116GHz and 
I c = 278.45 /iA and n = 1 ... 5. The agreement between 
the experimental data and this simple formula for the 
resonance condition is excellent. The effective capaci- 
tance of the junction calculated from the fitted zero bias 
plasma frequency and the critical current is found to be 
C = 1.61 pF. We note that, the corresponding classical 
harmonic resonances, arising from the nonlinearity of the 
junction, have been observed experimentally [p2| , p3| at 
large applied microwave power and at high temperatures. 

In Fig. [|, the escape rate T(I) reconstructed from the 
P(I) distribution is plotted for a range of microwave pow- 
ers for (a) the single- and (b) the two-photon absorption 
processes. In both cases, at low values of P MW the escape 
rate is a monotonic function of the bias current. With 
increasing P MW a clear resonance develops in the escape 
rate. The resonant current I r is indicated in the plot. 

The enhancement (r(P (tlw ) - r(0))/T(0) of the escape 



3 




FIG. 2: a) 3D plot and b) density plot of the measured P(I) distribution versus the applied microwave power P MW at v — 
36.554 GHz and T = 100 mK. c) Experimental data at 18.399 GHz for the same temperature. The switching probability P(I) 
is color coded as indicated by the scale. 



rate in the presence of microwaves of power P MW is plotted 
for both processes in the insets of Figs. ||a and b. P MW was 
chosen to result in a maximum enhancement of roughly 
10 in both cases. The data is fitted to a Lorentzian line 
shape H (solid line in Fig. |]). The width 8u of the single- 
photon resonance is in good agreement with the quality 
factor of the junction Q w v / Su — 380 determined from 
independent measurements |l6[ ] . For the two-photon pro- 
cess it is observed that the linewidth is approximately a 
factor of two larger than for the single photon process. 
For both processes it is observed that the line width of the 
transition is independent of the microwave power. This 
indicates that it is entirely limited by the life time of the 
excited state and that a coherent broadening of the res- 
onance was not observable in these measurements. The 
enhancement of the escape rate due to the microwave 
radiation increases approximately linearly with P MW for 
the single photon process, whereas it increases roughly 
as P/f w for the two-photon process. 
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FIG. 3: Applied microwave frequency v versus normalized 
resonant bias current I T /I C (dots). Dashed lines are a fit of 
the data to v = (1/n) u p (l - {I r /I c ) 2 ) 1/4 . 



We have compared our experimental results with the 
predictions of the Larkin-Ovchinnikov theory J24j . The 
bias-current dependent escape rate T(I) due to tunnel- 
ing from all possible energy levels was calculated using 
a master equation approach, considering the occupation 
of the energy levels in the presence of microwaves at 
T = 100 mK. The energy levels and matrix elements 
were determined using the approximations introduced 
in Refs. |^4|, p6| , which we found to be consistent 
with our direct numerical solutions of the Schrodinger 
equation for this problem. For single-photon processes 
the microwave-induced transition rates between nearest- 
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FIG. 4: Experimental escape rate V(I) for (a) single-photon 
and (b) two-photon absorption. Different curves correspond 
to P MW being increased (see arrow) from zero to a value at 
which the maximum enhancement F(P MW — r(0))/T(0) is ap- 
proximately 10. The resonance current I r is indicated by 
an arrow. The insets show the enhancement of the escape 
rate (r(P MW ) — r(0))/T(0) at the largest displayed value of 
P/iw Symbols are data, solid lines are fits to a Lorentzian 
line shape. 
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FIG. 5: P (I) distributions in the presence of microwave 
radiation at v = 36.554 GHz for P MW = — 7.5dB (open 
points) , 1 dB (gray points) and 4 dB (black points) measured 
at 100 mK. Dashed curves are calculated according to Larkin- 
Ovchinnikov theory. 



neighbor levels j are given by 
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where Tj is the inverse "lifetime" of the j + 1 — > J tran- 
sition. 

In Fig. § the measured switching current distributions 
for the single photon absorption are fitted to the calcu- 
lated distributions at different microwave powers for a 
junction capacitance of C = 1.27 pF, a critical current of 
I c = 278.25 /iA and an effective resistance of R = 500 51 
The capacitance obtained in this fit is about 20 % lower 
than that obtained from the plasma frequency. We sup- 
pose this discrepancy to be due to a microwave-induced 
level shift which is to be expected due to the large diag- 
onal matrix elements of the excited states. The theory 
well explains the power dependence of P(I), strength- 
ening the claim that the resonances are due to the mi- 
crowave induced transition of the phase from the ground 
state to the first excited state in the well. 

In the presented experiments we have found evidence 
for the microwave-induced multi-photon transitions be- 
tween quantized energy levels of the phase in a current- 
biased Josephson junction. This process could be yet an- 
other possible source of decoherence in microwave driven 
superconducting qubits. It could also be used for manip- 
ulating the quantum state of a qubit. 
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